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Abstract 

Until  recently  the  method  of  integral  operators  as  initiated 
by  S.  Bergman  [3]  and  I.  N.  Vekua  [75]  has  been  restricted  to  the 
case  of  elliptic  equations  and  the  investigation  of  steady  state 
phenomena.  In  these  lectures  we  survey  the  recent  developments 
on  the  use  of  integral  operators  to  investigate  equations  associated 
with  evolutionary  phenomena,  in  particular  parabolic  equations, 
pseudoparabolic  equations,  and  the  reduced  wave  equation  in  a 
stratified  medium.  The  topics  discussed  are  transformation 
operators  for  partial  differential  equations,  reflection  principles 
and  their  application,  the  propagation  of  radio  waves  around  the 
earth,  the  propagation  of  acoustic  waves  in  a spherically  stratified 
medium,  low  frequency  approximations  to  acoustic  scattering  problems 
in  a spherically  stratified  medium,  heat  conduction  in  two  temper- 
atures, inverse  problems  in  the  theory  of  heat  conduction,  and 
Runge's  theorem  for  parabolic  equations.  Open  problems  are  given 
at  the  end  of  each  section. 
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I.  Introduction 

Since  1960  a variety  of  books  ([3],  [5],  [10],  [11],  [41], 

[42],  [43],  [53],  [75],  [76]),  two  conference  proceedings  ([57], 

[66]),  plus  chapters  in  several  books  on  partial  differential 
equations  ([4],  [6],  [34],  [40],  [47],  [78])  have  been  written  on 

the  use  of  integral  operators  and/or  function  theoretic  methods 
in  the  theory  of  partial  differential  equations.  Most  of  this  work, 
until  recently,  has  been  concerned  with  elliptic  equations  and 
steady  state  problems.  Recently  an  analogous  approach  has  been 
discovered  for  equations  associated  with  evolutionary  phenomena, 
in  particular  parabolic  equations,  pseudoparabolic  equations,  and 
the  reduced  wave  equation  in  a stratified  medium.  These  lectures 
are  devoted  to  surveying  some  of  these  recent  developments,  and 
their  application.  Only  a brief  outline  of  proofs  will  be  given, 
and  the  reader  interested  in  detailed  proofs  is  referred  to  the 
list  of  references.  Although  most  of  the  material  presented  in 
these  lectures  has  been  developed  since  the  appearance  of  the  author's 
monographs  [10]  and  [11],  preliminary  versions  of  some  of  the  results 
of  these  lectures  can  be  found  in  these  books.  Open  problems,  of 
varying  difficulty,  will  be  given  at  the  end  of  each  section. 

Before  proceeding  perhaps  a biased  words  on  the  role  of 

integral  operators  in  the  theory  of  linear  partial  differential 
equations  are  in  older.  In  view  of  the  large  strides  that  have  been 
made  in  recent  years  in  developing  a general  theory  of  linear 
partial  differential  equations,  it  is  sometimes  suggested  that  the 
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theory  of  boundary  value  problems  for  linear  partial  differential 
equations  is  essentially  a closed  book.  That  this  is  far  from  the 
case  can  be  attested  to  by  any  applied  mathematician.  More  often 
than  not  the  general  theory  is  not  applicable  to  the  particular 
problem  being  investigated  due  to  the  fact  that  the  operator  is 
not  definite,  the  domain  is  unbounded,  the  boundary  data  is  dis- 
continuous, or  even  that  the  problem  is  improperly  posed  (c.f., 
Section  VTII) . In  other  cases,  problems  arise  due  to  the  need  to 
actually  compute  a solution  rather  than  "simply"  establish  its 
existence.  For  example,  an  existence  theorem  based  on  solving  an 
integral  equation  defined  over  an  unbounded  three  dimensional 
domain  is  often  of  limited  use  for  computational  purposes.  In 
order  to  handle  problems  such  as  these,  a variety  of  methods  have 
been  developed  which  are  roughly  speaking  characterized  by  the  fact 
that,  although  they  are  highly  effective  for  the  problems  they  are 
designed  to  treat,  by  their  very  nature  are  restricted  to  rather 
limited  classes  of  equations.  Perhaps  the  best  example  of  this  is 
the  use  of  generalized  double  and  single  layer  potentials  in  the 
study  of  partial  differential  equations  with  constant  coefficients 
(although  this  method  can  be  applied  to  equations  with  variable  co- 
efficients, the  practicality  of  such  an  approach  from  the  point  of 
view  of  analytic  and  numerical  approximations  is  rather  limited) . 
From  one  point  of  view  the  use  of  generalized  potential  theory  can 
be  seen  as  a branch  of  the  theory  of  integral  operators,  i.e.,  those 
operators  mapping  continuous  functions  onto  solutions  of  linear 


partial  differential  equations  with  constant  coefficients.  However, 
in  these  lectures  integral  operators  will,  in  general,  be  viewed  in 
a more  restrictive  sense,  i.e.,  that  in  which  the  integral  operator 
maps  solutions  of  partial  differential  equation  with  constant 
coefficients  onto  solutions  of  a partial  differential  equation  with 
variable  coefficients  (or  occasionally  analytic  functions  onto 
solutions  of  a partial  differential  equation) . Obviously,  integral 
operators  in  the  latter  sense  and  generalized  potential  theory  can 
often  be  combined  (c.f..  Section  V),  and  hence  the  distinction  is 
not  a sharp  one.  In  fact  it  is  perhaps  more  fruitful  to  view 
generalized  potential  theory  and  the  method  of  integral  operators 
as  complimentary,  one  dealing  with  partial  differential  equations 
with  constant  coefficients  and  the  other  treating  certain  classes 
of  partial  differential  equations  with  variable  coefficients.  In 
view  of  the  above  described  role  of  the  method  of  integral  operators 
in  the  theory  of  linear  partial  differential  equations,  the  primary 
interest  is  not  so  much  in  the  integral  operator  itself,  but  rather 
in  how  it  can  be  used  to  yield  constructive  methods  for  solving 
boundary  value  problems  arising  in  various  areas  of  applied  mathe- 
matics. Of  course,  as  with  any  fruitful  area  of  mathematics,  in 
the  process  of  achieving  the  desired  aim  many  other  problems  of 
independent  interest  arise  along  the  way.  Hence,  although  our  pri- 
mary aim  in  these  lectures  is  to  illustrate  the  use  of  integral 
operators  in  applied  mathematics  by  considering  various  problems 
arising  in  the  theory  of  wave  propagation  and  heat  conduction,  we 
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shall  also  find  the  need  to  examine  such  topics  as  reflection  prin- 
ciples and  the  analytic  continuation  of  solutions  to  partial  dif- 
ferential equations,  Range's  theorem  for  parabolic  equations, 
generalized  moment  problems,  and  the  completeness  of  certain  systems 
of  entire  functions. 


II.  Transformation  Operators  for  Partial  Differential  Equations  . 

In  this  section  we  shall  construct  the  integral  operators  to 
be  used  in  the  next  four  sections.  The  term  transformation  opera- 
tors is  used  because  of  the  similarity  of  our  operators  to  trans- 
formation operators  for  ordinary  differential  equations  (c.f., 
Appendix  4 of  [56]).  We  first  consider  the  parabolic  equation 


u + a(x,t)u  + b(x,t)u  = c(x,t)u,  ; c(x,t)>  o (2.1) 

and  note  that  by  a change  of  independent  and  dependent  variables 
(2.1)  can  be  reduced  to  an  equation  of  the  form 

u + q(x,t)u  = u,  . (2.2) 

XX  ^ t 


Hence  we  restrict  our  attention  to  (2.2)  and  for  the  salce  of  sim- 
plicity assume  that  q(x,t)  is  an  entire  function  of  x and  t . 
We  loo)<  for  solutions  of  (2.2)  defined  in  a neighborliood  of  the 
origin  in  the  form 


u ( X , t ) = li  ( X , t ) + 


fX 


-^x 


P(s,x,t)  h(s,t)ds 


(2.3) 


(I  t P)h 


5. 


where  h(x,t)  is  a solution  of  the  heat  equation 

h = h 

XX  t 

Substituting  (2.3)  into  (2.2)  and  integrating  by  parts  show 
that  (2.3)  will  be  a solution  of  (2.2)  if 


P - P 
XX  ss 

P(x,x,t)  = 


+ q (x, t) P = P 


rx 


q (s,t) ds 


P(-x,x,t)  = 0 


(2.5) 


A solution  of  (2.5)  can  be  constructed  by  iteration  (c.f.  [11],  [12]). 

Note  that  since  P is  a Volterra  operator,  1 + P is  invertible. 

For  future  applications  we  note  that  if  u(0,t)  = 0 then  (2.3)  can 
be  rewritten  as 


u (x, t) 


h (x  , t)  + 


X 

[P(s,x,t) 

0 


P(-s,x,t) ]h(s, t)ds 


= h (x, t)  + 


rX 

P(s,x,t)h(s,t)ds. 

0 


(2.6) 


By  using  the  operator  I + P many  properties  of  the  solutions  to  the 
heat  equation  can  also  be  obtained  for  solutions  of  (2.2).  As  an 
example  we  have  the  following  theorem: 

Theorem  2.1:  Let  be  a sequence  of  complex  numbers  such  that 


1 im 

n-*-oo 


> 


0 


Then  {(I+P)h~}  is  complete  (with  respect  to  the  maximum  norm) 
in  any  rectangle  R for  solutions  u(x,t)  of  (2.2). 

Proof : Let  u(x,t)  be  a solution  of  (2.2)  in  R such  that  u(x,t) 

is  continuous  in  R . We  can  write  u = (I+P)h.  Hence  it  suffices 
to  show  that  {h  } is  a complete  family  of  solutions  for  h = h. 
in  R.  In  R,  h(x,t)  can  be  approximated  by  an  entire  solution 
h(x,t)  of  the  heat  equation  ([11],  [15]),  and  h(x,t)  depends 


continuously  on  its  Cauchy  data  in  any  thin  complex  neighborhood  of 

/-  2 ^ 

■ 

X = 0 ([11])  . The  result  now  follows  from  the  fact  that  le  / 


are  complete  in  such  a neighborhood  if  lim  — ^ > 0 ([56]). 

n - 

We  now  turn  our  attention  to  elliptic  equations  with  spherically 
symmetric  coefficients  defined  in  a domain  in  Euclidean  three  space 

3 

R . We  first  consider  elliptic  equations  defined  in  bounded  domains 
and  present  a version  of  R.  P.  Gilbert's  "method  of  ascent"  ([42], 
[44]).  The  construction  given  below  is  due  to  D.  Colton  and  R.  Kress 


([31]).  Let  D be  a bounded  starlike  domain  in 
D C.  {(x|  < a}  and  consider  the  equation. 


such  that 


A^u  + X (1+B(r) ) u = 0 


(2.7) 


defined  in  D where  B(r)  e C [0,a],  x = ]x|.  We  look  for  a 


solution  of  (2.7)  in  D in  the  form 


.7 


u(r,0,i))) 


h ( r , 0 , ij) ) 


+ 


r 

G{r,s;A)h(s,9,4))ds 

0 


= (I+G)h 

where  (r,0,4))  denote  spherical  coordinates  and 


(2.8) 


A^h  = 0 


(2.9) 


in  D.  Substituting  (2.8)  into  (2.7)  and  integrating  by  parts 
shows  that  (2.8)  will  be  a solution  of  (2.7)  provided 


r^[G 


rr 


+ I + A^H-B(r))Gl  = + | GJ 


G(r,r; A) 


r 

s (1+B (s) ) ds 

0 


(2.10) 


and  G(r,0;A)  is  bounded  for  0 £ r _<  a.  The  solution  of  (2.10)  can 
be  found  by  iteration  ([31]).  The  function  G(r,s;A)  can  in  fact 
be  related  to  the  Piemann  function  for  a related  hyperbolic  equation; 

c 1/2 

Theorem  2.2:  G(r,s;A)  = - (p)  R^(r,r;s,0)  where  R(x,y;J;,n)  is 

A 2 

the  Rieinann  function  for  R + — j (l+B(/x^))R  = 0 where  the 

xy  ^ 

subscript  denotes  differentiation  with  respect  to  C- 


Proof : [42],  [44]. 

We  now  want  to  present  a complimentary  operator  to  I + G valid 
for  exterior  domains.  Tliis  operator  is  due  to  D.  Colton  and 
W.  Wendland  ((33))  and  0.  Colton  and  R.  Kress  ([30]).  Let  D be 


. 8 


a bounded  star like  domain  in  E and  b a real  number  such  that 


{lx|  < b}  d D.  Consider  (2.7)  defined  in  E’^\D  and  assume  that 


2 

B(r)  e C^(b,“),  B(r)  = 0{e  ^ ),  and  A < 2yh.  We  look  for  a 


3. 


solution  of  (2.7)  in  E' \d  in  the  form 


u(r,0,4))  = h(r,0,iji)  + 


K(r,s;A)h(s,0,()))ds 


(2.11) 


= (I+K)h 


where  now  h(r,0,(t))  is  a solution  of 


A.jh  + A h = 0 


(2.12) 


.3^ 


in  E \d.  It  is  also  possible  to  construct  an  operator  whose 
do.main  is  solutions  of  Laplace's  equation;  however  for  purposes  of 
application  to  problems  in  scattering  theory  we  restrict  ourselves 
here  to  an  operator  whose  domain  is  solutions  of  the  Helmholtz 
equation.  Substituting  (2.11)  into  (2.7)  and  integrating  by  parts 
shows  that  (2.11)  will  be  a solution  of  (2.7)  provided 


r^[K  + - K + A'"  (H-B(r)  )K]  = s^  [K  + - K + A^K] 
r r r r s s s s 


K ( r , r ; A ) = - 


2r 


SB ( s ) ds 


(2.13) 


K(r,s;A)  - 0(e  YJ-'stA  s/4Yr^  b £ r £ 


The  solution  of  (2.13)  can  be  found  by  iteration  ((30]). 


> 


so  does  u(r,0,(())  = (r+K)h.  We  finally  observe  that  since  G and 
are  Volterra  operators,  I + G and  I + K are  invertible. 


Open  Problem;  Obtain  asymptotic  estimates  as  A -*■  +»  for 
G{r,s;A)  and  K(r,s;A).  Such  results  would  be  useful  in  connection 
with  various  problems  arising  in  scattering  theory. 


Open  Problem;  Construct  transformation  operators  for 

A^u  + B(r)u  = u^  having  as  domains  solutions  of  the  heat  equation 

A^u  = u^.  (For  partial  progress  in  this  direction  see  [65]). 


111.  Reflection  Principles  and  Their  Applications. 


We  first  consider  the  parabolic  equation 


u , + a(x,t)u  + b(x,t)u  = u 

XX  X u 


(3.1) 


defined  in  a domain  of  the  form  pictured  below. 


.10 

x=S2 (t) 


2 0 — 

We  assume  that  u(x,t)  e C (D)  C (D)  , that  a(x,t)  and  b(x,t) 
are  entire  functions  of  x and  t (this  assumption  can  be  relaxed) , 
and  that  s^(t)  and  ^2^^^  analytic  in  a neighborhood  of 

(0,tQ].  Define  the  "reflection"  of  D across  a : x = s^(t)  by 

D*  = {(x,t)  : 2s^(t)-S2(t)  < x < s^(t),  0 < t < tp}. 

Theorer-  3.1;  Suppose  u(x,t)  is  a solution  of  (3.1)  such  that 
u(x,t)  = f(t)  on  a where  f(t)  is  analytic  in  a neighborhood 
of  [Ojtg].  Then  u(x,t)  can  be  uniquely  continued  as  a solution 
of  (3.1)  into  D (J  a U D . 

Proof  ([11],  [12],  [13]);  A change  of  variables  transforms  the 

problem  into  the  case  when  a(x,t)  = 0,  Sj^(t)  =0.  By  constructing 
a special  solution  of  a non-characteristic  Cauchy  problem  for  (3.1) 
we  can  further  reduce  the  problem  to  the  case  when  f't)  = 0.  From 
Section  II  we  have 

•X  ^ 

u(x,t)  = )i(x,t)  + P (s  , X,  t)  h (s  , t)  ds 

0 

where  h(0,t)  = 0.  The  Theorem  now  follows  from  the  reflection 


'i 


principle  Cor  the  licat  equation  ([80]). 
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Theorem  3.2:  Theorem  3.1  remains  valid  if  the  condition 
u(x,t)  = f{t)  on  a is  replaced  by 

a(t)u(x,t)  + 6(t)u^(x,t)  + Y(t)u^(x,t)  = f(t)  on  a 

provided  a(t),  B(t),  Y(t)  and  f(t)  are  analytic  in  a neigh- 

borhood of  [Ojtp]  and  ti(t)  = (B(t),Y(t))  is  never  tangent  to 
0 and  either  never  parallel  to  the  x-axis  or  always  parallel  to 
the  x-axis. 


Proof ; [14 ] . 

The  above  reflection  principles  can  be  used  to  extend  Theorem 
2.1  from  the  case  of  solutions  of  (2.2)  defined  in  a rectangle  to 
the  case  of  solutions  of  (2.2)  defined  in  a domain  with  moving 
boundaries : 


Theorem  3.3:  The  set  {(I+P)hj^}  is  a complete  family  of  solutions 
to  (2.2)  in  D for 

^ ±iA  x-A^t 

1)  h“(x,t)  = e " , lim  > 0, 

" n-voo  X 2 

n 


or 


2)  hj^(x.t) 


[n/2]  n-2k^k 
Y ^ t 

kio  <'^-^X)!k!  • 


Proof  ( [ 11 1 , [16  1);  Approximate  the  boundary  data  by  analytic 

functions,  reflect  repeatedly  across  x = Sj^(t)  and  x = s^(t), 
and  use  the  fact  that  the  sets  1)  and  2)  defined  above  are  complete 
for  solutions  of  the  heat  equation  defined  in  a rectangle. 


.12 


For  numerical  applications  of  Theorem  3.3  see  [8]. 

We  now  turn  our  attention  to  deriving  a reflection  principle 
for  solutions  to  the  Helmholtz  equation  in  which  is  analogous 

to  the  Schwarz  reflection  principle  for  harmonic  functions 
vanishing  on  a portion  of  a spherical  boundary.  Our  aim  is  to 
then  use  this  reflection  principle  to  deduce  a continuation 
theorem  connected  with  the  inverse  scattering  problem  for  acoustic 
waves.  Let  D be  a bounded,  starlike  domain  in  containing 

the  ball  S = {x  : |xl  £ a}  in  its  interior,  and  let 
u(x)  = u(r,0),  0 = (0  , . . . , , be  a solution  of 

A u + A^u  = 0 (3.2) 

n 

in  D\S . 


If  u(a,0)  = 0,  then  in  D\S  we  can  represent  u(r,0)  in  the 
form  ([27]) 


u(r,0)  = h(r,0)  + 


s^  (r,s; A)h (s, 0) ds 


(3.3) 


f 


I 

f 


A h = 0 
n 


(3.4) 


in  D\S,  h(a,0)  = 0,  and  K(r,s;X)  is  a solution  of  the  initial 
value  problem 


K ( r , r ; X ) 


2-n,  2 2. 

r (r  -a  ) 


(3.5) 


K(r,a;X)  = 0 . 

The  solution  of  (3.5)  can  be  found  by  iteration  ([27]).  Hence 
we  have  the  following  theorem; 

Theorem  3.4;  If  u(r,0)  is-a  solution  of  (2.2)  in  D\S  such 

that  u(a,0)  = 0,  then  u(r,0)  can  be  uniquely  continued  as  a 

★ 

solution  of  (2.2)  into  D = 

Proof  ( [27] ) ; This  follows  from  the  representation  (3.3)  and  the 
Schwarz  reflection  principle  for  harmonic  functions. 

Now  let  n = 3 and  suppose  u(x)  = u(r,0,(ti)  is  a solution 
of  (3.2)  in  the  exterior  of  a bounded  simply  connected  domain 
n such  that  u(r,0,4')  satisfies  the  Sommerfeld  radiation  condition 

lim  r (Uj,-iXu)  = 0 . 

f-VOO 


(3.6) 
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d 


Then  at  infinity 

u ( r , 0 , ij) ) 


iXr 

^-7-  f(e,4>;X) 


(3.7) 


i X IT  c o 5 0 

Suppose  that  on  we  have  u(r,0,0)  = e . The  inverse 

scattering  problem  is  to  determine  Q from  a Icnowledge  of 
f(0,<t);X).  It  is  known  that  f(0,(J>;X)  uniquely  determines  n 
([55]).  From  f(0,d);X)  we  can  determine  u(r,0,4))  outside 
the  smallest  ball  S containing  in  its  interior  (c.f.  [60]). 

Hence  the  problem  is  to  continue  u(r,0,(}))  across  9S  and 
determine  the  locus  of  points  such  that  u(r,0,(J))  - = 0. 

This  can  be  done  (in  theory!)  by  using  Theorem  3.4  (after  a 
reduction  to  the  case  when  u(r,0,(t))  =0  on  3S  - c.f.  [27]) 
in  conjunction  with  R.  P.  Gilbert's  "envelope  method"  (c.f.  [10], 

[41],  [42])  as  applied  to  harmonic  functions. 

Definition  3.1;  Let  f(z)  be  an  entire  function  of  exponential 
type.  Then  the  indicator  diagram  of  f(z)  is  the  interior  of 
the  convex  hull  of  the  singularities  of  the  Laplace  transform 
of  f(z)  (The  indicator  diagram  can  also  be  defined  in  terms 
of  the  growth  of  f(z)  along  rays  passing  through  the  origin  - 
c.f.  [56]). 

Theorem  3.5;  Let  u(r,0,(())  = u(r,0)  be  axially  symmetric  (which 
implies  f(0,()))  = h(cos  0)).  Define 


F(z) 


(l+4z^)d/: 

(l-4iz^-4z^7^^ 


(Note  that  h(0 
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CO  oo 

= I a^P^(0  F(z)  = 2 I a^(2iz)").  Then 
n=0  ^ n=0  " 

F(z)  can  be  continued  to  an  entire  function  of  exponential  type. 
If  I is  the  indicator  diagram  of  F(z),  then  u(r,6)  is 
regular  in  the  exterior  of  lU  I,  where  the  bar  denotes  complex 
conjugation . 

Proof;  [101,  [24]. 

Remark ; In  [lOJ  and  [24]  it  was  not  possible  to  exclude  the 
possibility  that  u(r,0)  was  singular  along  the  axis  9=0  or 
6 = TT  since  at  that  time  Theorem  3.4  was  unavailable.  However, 
if  in  the  proofs  of  [10]  and  [24]  one  uses  Theorem  3.4  instead 
of  Lewy's  reflection  principle  as  applied  to  the  axially  symmetric 
Helmholtz  equation,  one  arrives  at  Theorem  3.5. 

For  another  approach  to  the  inverse  scattering  problem  see 
[79]. 

Open  Problem;  Derive  an  analogue  of  Theorem  3.5  when  Q is  no 
longer  a bounded  domain. 

Open  Problem;  In  Theorem  3.2  remove  the  restriction  that  ij(t] 
is  either  never  parallel  to  the  x-axis  or  always  parallel  to  the 
x-axis . 

IV.  The  Propagation  of  Radio  Waves  Around  the  Earth. 

In  this  section  we  shall  show  how  the  transformation 
operators  obtained  in  Section  II  can  be  used  to  construct  approxi- 
mate solutions  to  the  Fock-Leontovich  equations  describing 
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radiowave  propayation  around  the  earth  under  the  assumption  of 
a spherically  stratified  atmosphere  but  ignoring  the  effect  of 
the  ionosphere  ([38]).  We  assume  that  a vertical  electric  dipole 
is  situated  on  the  surface  of  the  earth,  and  without  loss  of 
generality  we  assume  that  this  point  is  the  north  pole.  Let 
(x,y)  denote  the  point  of  observation  of  the  electromagnetic 
field,  where  x is  the  (normalized)  distance  along  the  surface  of 
the  earth  from  the  north  pole  to  a point  directly  below  the 
observation  point,  and  y is  the  (normalized)  distance  of  the 
observation  point  to  the  surface  of  the  earth.  Let  a denote 
the  radius  of  the  earth,  k the  wave  number,  and  w(x,y)  the 
(normalized)  Hertz  potential  (For  precise  definitions  the  reader 
is  referred  to  [38]).  Then  under  the  assumption  of  a perfectly 
conductive  earth  and  using  the  fact  that  ka  is  very  large,  we 
are  led  to  the  following  initial-boundary  value  problem  for 
w(x,y)  : 

Wyy  + iw^  + y(l+g(y))w  = 0 (4.1a) 


Wy(x,0)  = 0;  0 < X < Xq 


(4.1b) 


2 2 

= ~ expdy^  /4x)  ; 0 < x < x^ 


(4.1c) 


1 im 

X '0 


w(x,y> 


— ^ exp (iy^/4x) 
/x 


= 0; 


0 < y < yg 


(4. Id) 


where  g(y)  is  a real-valued  slowly  varying  function  (in  our 
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case  assumed  to  be  a polynomial)  related  to  the  refractive  index 
of  the  atmosphere  and  are  positive  constants  such  that 

X < Xq , y < y^ , y^  >>  (4.1a)  is  the  Fock-Leontovich 

equation,  (4. lb)  reflects  the  fact  that  the  earth  is  perfectly 
conductive,  (4.1c)  is  a consequence  of  the  Fresnel  reflection 
law  and  geometrical  optics,  and  (4. Id)  is  due  to  the  presence  of 
a vertical  electric  dipole  at  the  north  pole.  For  a full 
discussion  of  the  derivation  of  ( 4 . la) - ( 4 . Id)  see  [38]  and  [17]. 
The  advantage  of  studying  the  equations  (4 . la) - (4 . Id)  instead  of 
the  usual  boundary  value  problem  associated  with  the  Helmholtz 

equation  is  that  it  permits  the  investigation  of  not  only  the 

2 2 
illuminated  region  (y  >>  x ) and  the  shadow  region  (y  <<  x ) 

but  also  intermediate  cases,  namely  the  region  of  the  penumbra. 


vertical  dipole 


shadow 

region 


illuminated  region 


penumbra 


iliK 


earth 


We  look  for  a solution  of  ( 4 . la) - ( 4 . Id)  in  the  region 
D = { (x,y)  : 0 < X < Xq , 0 < y < y^}  in  the  form 


L 


2 2 

w(x,y)  = — exp(iy  /4x) 


+ — [P (s , y) +P (-S ,y) ] exp  (is  /4x)ds  + u(x,y) 
/x  0 


(4.2) 


where  P(s,y)  is  the  kernel  P(s,y,t)  defined  in  Section  II 

2 

(which  in  our  case  is  independent  of  t)  and  u(x,y)  e C (D) 
/^C^(D)  is  a solution  of  (4.1a)  such  that  Uy(x,0)  = 0.  We 
now  want  to  determine  the  remaining  initial-boundary  data  for 
u(x,y)  and  derive  a method  for  approximating  u(x,y)  as  well 
as  the  term 


2 fy  2 

— (P (s , y) +p (-S ,y) ] exp ( is  /4x)ds 

Jo 


(4.3) 


appearing  in  (4.2).  Note  that  this  last  approximation  is  non- 
trivial due  to  the  factor  of  2//x  appearing  in  front  of  the 
integral  sign.  In  view  of  the  various  approximations  which  need 
to  be  made  in  order  to  determine  the  initial-boundary  data  for 
u(x,y),  as  well  as  the  fact  that  the  boundary  data  (4.1c)  is  an 
approximation  to  begin  with,  the  following  a priori  estimate  is 
of  basic  importance; 

Theorem  4.1;  Let  u(x,y)  f-  C (D)  C (D)  be  a solution  of  (4.1a) 
in  D such  that  Uy(x,0)  = 0,  u(x,yQ)  = f(x),  u(0,y)  = h(y). 
Then  there  exists  a positive  constant  M such  that 


fs 


s 2 f 0 2 r ® 2 

u (t , n) d r dsdn  £ M ) f (s)  | ds  + ) h (n)  | dn 

0 •'0  ■’  0 


0 


0 
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Proof  ([17]) ; Define 


(x,y) 


X 

u (t ,y) dT 

0 


X 

f (t ) dt . 

0 


Then  u^{x,y)  satisfies  a non-homogeneous  version  of  (4.1a) 

with  homogeneous  initial-boundary  data.  The  result  now  follows 

by  applying  standard  eigenfunction  expansion  methods.  Note  that 

2 

due  to  the  fact  that  in  general  f'(x)  / L (0,Xq),  it  is  not 

2 

possible  (by  these  methods)  to  replace  the  weighted  L -norm 
2 

by  a L norm  on  the  left  hand  side  of  the  inequality  in  Theorem 

4.1. 


We  now  proceed  to  approximate  (4.3)  and  the  initial-boundary 
data  f(x)  and  h(y)  as  defined  in  Theorem  4,1.  From  a result 
on  asymptotic  expansions  due  to  Erdelyi  ([36])  it  is  possible  to 
show  that 


_2 


2 

[P (s,y) +P {-s,y) ] exp (is  /4x)ds 


2/2TT(l+i)P(0,y)+R(x,y) 

(4.4) 


where  |R(x,y)|  < constant • max  |P  (s , y) -P  (-s ,y ) | . By  using 

D ® ® 

this  result  and  the  asymptotic  properties  of  Fresnel  integrals 

(after  approximating  P (s ,y) +P  (-s ,y)  by  a polynomial)  one  can 

2 1/2 

shov/  that  to  within  0(4xQ/yQ  ) we  have 
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f (x)  = Pj^(x)  + exp  ( iy‘‘/4x)  P2  (x) 
h(y)  = P3(y) 


(4.5) 


.19a 


where  (x)  , ^2^^^  polynomials.  A similar 

analysis  allows  us  to  approximate  (4.3)  for  0 £ x £ Xq , 

0 £ y £ y^ . For  details  the  reader  is  referred  to  [17]. 


Our  task  now  is  to  approximate  u(x,y)  satisfying  (4.1a) 
along  with  the  initial-boundary  data  Uy(x,0)  = 0,  u(x,yQ)  = f (x) , 
u(0,y)  = h(y)  where  f(x)  and  h(y)  are  given  by  (4.5).  To 
this  end  we  shall  again  make  use  of  the  transformation  operators 
constructed  in  Section  II. 


2 

Theorem  4.2:  Assume  x_  < y_  and  let  A 
0 — 0 n 

Let 


2mr 


(n-1/2)^. 

^0 


-iAnX 

h2n(^'y)  = cos/A^  ye  , n = ±1,±2,... 

-ip  ^x 

h2n+i(^'Y)  = cos  M^y  e , n = 1,2,... 


Then  with  respect  to  the  norm  defined  in  Theorem  4.1,  the  set 
{Uj^}  where 


u^(x,y) 


= hj^(x,y)  + 


y 

[P(s,y)+P(-s,y) ]h(s,y)ds 

0 


is  a complete  set  of  solutions  to  (4.1a)  satisfying  (4.1b). 


Proof  ([17]): 

set  (h  } is 
n 

hy(x,0)  = 0. 


As  in  Theorem  2.1,  it  sufficies  to  show  that  the 
complete  for  solutions  of  h^y+ih^  = 0 satisfying 
But  this  follows  from  Theorem  4.1  and  Levinson's 
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f -iA^xj  2 

result  that  \e  J is  complete  in  L (0,Xq)  if 

Uni  1 ~ (Ul  +1)  (c.f.  [56]). 

To  approximate  u(x,y)  we  now  set 

N 

u^(x,y)  = I aj^Uj^(x,y) 


(4.6) 


. 2N 

for  a given  integer  N and  minimize  (in  I ) the  quadratic 

functional 


Q(a  ,...,a  ) = I I u^'^  (x,y  ) -f  (x)  I U, 

^ M u IIl‘-(0,Xq) 

(4.7) 

+ I |u^(0,y)-h(y)  I U- 

" 11^(0, y^)  . 


Open  Problem:  Treat  the  case  when  the  earth  is  no  longer  perfectly 
conductive.  (The  generalization  is  non-trivial ! ) 

Open  Problem:  Treat  the  case  when  i’Q  = “»  i.e.  avoid  the  use 
of  the  geometric  optics  approximation  (4.1c). 

V . The  Propagation  of  Acoustic  Waves  in  a Splierically 
S bra ti lied  Medium. 

In  this  section  we  shall  show  how  the  transformation  operators' 
constructed  in  Section  II  can  be  used  in  the  investigation  of 
acoustic  scattering  problems  in  a spherically  stratified  medium. 

We  first  consider  the  case  when  the  incoming  wave  is  scattered  by 


1 


j 


1 

! 


I 


the  presence  of  a spherically  stratified  quasi-hcmogeneous 
medium,  but  in  the  absence  of  any  obstacle.  If  we  denote  the 
velocity  potential  by  u(x^  and  the  velocity  potential  of  the 
scattered  wave  by  (factoring  out  a term  of  the  form 

^-loitj  arrive  at  the  following  set  of  equations  for  the 

determination  of  u^lx): 


u{x)  = + Ug(x) 


A^u  + X^(l+B(r))u  =0  in 


(5.1) 


lim  ^ ~ 


where  A is  the  wave  number,  x = (x,y,z),  r = |x|,  and 


B(r)  - 


^0  P 


- 1 where  c(r)  is  the  speed  of  sound  and 


— Yir 

c-  = lim  c(r).  We  make  the  assumption  that  B(r)  = 0(e  ) 

r-f-oo 

and  0 < A < 2y(/2-1) . We  look  for  a solution  of  (5.1)  in  the 
following  form:  For  r < 1 we  represent  u(x)  as 


u(x)  = ) a u (r)P  (cos  0 

In  >■' 
n=0 


u^^(r)  = (I+G)r 


(5.2) 


where  I + G is  the  transformation  operator  defined  in  Section  II, 
P^(cos  0)  is  Legendre's  polynomial,  and  the  constants  a^  are  to 
be  determined.  For  r > 1 we  represent  u(x)  as 


be  determined. 
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u (x) 


(I+K) 


+ y b (Xr)P  (cos 

n=0  ^ 


(5.3) 


where  (I+K)  is  the  transformation  operator  defined  in  Section 

II,  (Xr)  is  a spherical  Hankel  function,  and  the  constants 

b are  to  be  determined.  If  we  now  use  Sonine's  formula  to 

n 

expand  e^^^  and  require  u(x)  and  its  first  derivatives  to 
agree  at  r = 1 we  are  led  to  the  following  algebraic  system  for 
the  determination  of  the  constants  a^^  and  b^: 


a u (1)  - b 
n n n 


(I+K)hj^^^  (Xr) 


= (2n+l)i 


r=l 


(I+K) j^(Xr) 


r=l 


a u ' (1)  - b^ 
n n n 


^ (I+K)h  (Xr) 
dr  ~ ~ n 


= (2n+l)i 


r=l 


dF 


(5.4) 


r=l 


where  j^(Xr)  denotes  a spherical  Bessel  function.  From  (5.4) 
it  is  now  possible  to  deduce. by  using  Crammer's  rule  and  uniform 
estimates  for  Bessel  and  F^an)cel  functions  that 


K(l,s;X)j^(Xs)ds 


O 


■fx 

n- 

[l2gj 

K(l,s;X)h^^^  (Xs)ds|  = 0 


r(n+l/2)2'^ 

nX 


(5.5) 


where  g - y - X^/4y  and  K(r,s;X)  is  the  l^ernel  of  the  operator 
I+K.  These  estimates  imply  that 


which  imply  that  the  series  representations  for  u{x)  converge 
and  satisfy  (5.1).  The  uniqueness  of  the  solution  follows  from 
the  fact  that  for  X real,  (I+K)h^^^(Xr)  are 

linearly  independent  solutions  of 

X^(l  + B(r))-  y = 0.  (5.7) 

r 

For  details  of  the  above  calculations  the  reader  is  referred  to 
( 26 ] and  [ 30 ] . 

We  now  turn  our  attention  to  the  same  problem  considered 
above,  except  that  in  addition  to  the  spherically  stratified 
medium  there  is  a "hard"  obstacle  D present.  We  assume  that  D 
is  bounded,  starlilce,  and  has  smooth  boundary  90  with  outward 
pointing  unit  normal  v.  In  this  case  the  equations  for  deter- 
mining the  velocity  potential  u(x)  become 


Du  . 

lim  r 

Dr 


0 


(5.8) 
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We  again  assume  that  B(r)  = 0(e  ) and  restrict  ourselves  to 
values  of  X such  that  X < 2Yb  where  {|x|  < b}CZ  D.  We  look 
for  a solution  in  the  form 


u(r,0,i|i)  = h(r,G,(^)  + K ( r , s ; X)  h ( s , 0 , (Jj  ) ds 


(5.9) 


where  K(r,s;X)  is  the  kernel  of  the  transformation  operator 
I + K and  h(x)  is  a solution  of 


A^h  + X h = 0 


(5.10) 


in  E \D  of  the  form 


h (x)  = e^^^  + hg (x) 


lim  r hr:—  - iXh  = 0 


(5.11) 


Following  Ursell  ([74];  see  also  [50])  we  represent  in 

the  form 


h^(x)  = ij(C)G(5;x;X)da)^ 


(5.12) 


where  m(C)  is  a potential  to  be  determined  and  G(f;;x;X)  is 
a fundamental  solution  of  the  Helmholtz  equation  in  the  exterior 
of  {[xl  £ Ij]  satisfying  the  Sommerfeld  radiation  condition,  and 
on  I X 1 = b 


-r  - 4 - a G = 0 
f>  r 


(5.13) 


where  a 


0 < 6 < TT.  The  function 


G(C;x;A)  can  be 
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i6 
e , 


constructed  by  separation  of  variables  and  is  introduced  to  avoid 

the  problem  of  the  non-invertibility  of  the  integral  equation 

associated  with  scattering  problems  for  X an  eigenvalue  of  the 

interior  Dirichlet  problem  of  the  Helmholtz  equation  ([74]).  if 

iXz 


+ h,(x)  into  (5.9)  with  h (x)  defined 
s ~ s ~ 


we  now  substitute  e 
by  (5.12),  interchange  orders  of  integration  and  let  x tend  to 
9D,  we  arrive  at  a Fredholm  integral  equation  for  the  determination 
of  u(^)  of  the  form 


^ f(x)  = (I+T(X))y  ; X e 9D  (5.14) 


where  f (x)  = 

(I+K)e^^^ 

and  T(X)  is  a compact 

integral  operator 

on  the  Banach 

space  OD) 

([11],  [33]). 

Theorem  5.1: 

Let  0 < X < 

2Yb.  Then  in  C*^OD), 

(I+T(X) )"^ 

exists . 

Proof  ( [3  3] ) : 

Since  K is 

a Volterra  operator,  (I+K)  ^ exists. 

and  hence  if 

(I+T(X))ii  = 0 

we  can  conclude  that 

(I+To(X)  )y  = 0 

vfhere  Tq  = T 

From 

the  choice  of  G(5;x;X) 

we  can  now 

II 

o 

conclude  that 

11  ( O = 0,  arid  hence  by  the  Fredholm 

Alternative 

( I+T (X) ) ^ exists . 

[Jote  tiiat  one  advantage  of  the  integral  operator  method  outlined 
above  over  the  straiglit  forward  use  of  integral  equations  is  that 
the  resulting  integral  equation  is  defined  over  a bounded  two 
dimensional  region  instead  of  an  unbounded  three  dimensional 

^rjegion.  


Open  Problem;  Investigate  problem  (5.8)  in  the  case  when  D 
is  no  longer  starlike  with  respect  to  the  origin.  For  a pre- 
liminary investigation  of  this  problem  see  [29] . 
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Open  Problem;  Viewing  I!T(X)  as  an  operator  valued  analytic 
function  of  A,  determine  the  nearest  singularities  to  the 
origin  of  (IfT(A))  Such  a result  would  lead  to  new  construc- 

tive methods  fc-r  solving  15,8)  obtained  by  expanding  (I+T(A))  ^ 
in  powers  of  A and  determining  the  coefficients  recursively. 

For  a survey  of  the  results  of  this  section  in  the  case 
when  B(r)  has  compact  support  see  [28], 


VI.  Low  Frequency  Approximations  to  Acoustic  Scattering 
Problems  in  a Spherically  Stratified  Medium. 


In  this  section  we  shall  consider  the  problem  of  approximating 
the  solutions  of  the  problems  discussed  in  the  previous  section 
in  the  case  when  the  wave  number  A is  small.  In  the  case  of  a 
homogeneous  medium,  the  problem  of  low  frequency  approximations  to 
scattering  problems  has  been  extensively  investigated  by  Ralph 
Kleinman  and  his  co-workers  (c.f.  [51],  152]).  We  first 

consider  the  problem  of  scattering  by  a spherically  stratified 
medium  in  the  absence  of  an  obstacle,  i.e.  to  determine  the  velocity 
potential  u(x)  from  the  equation 


u ( X ) = e ^ ^ s ^ ^ ^ 


A^u  + A (1+B(r))u  = 0 in 


<lu, 

aT 


- - i A u 


] im  r 


0 


(6.1) 
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where  we  now  assume  that  X,  in  some  sense,  is  small.  In  this 
case  a classical  method  for  solving  (6.1)  is  to  reformulate 
the  problem  as  the  integral  equation 


I \ 1 2 

u{x)  = e + X 


J 

.3 

4it  |x-^| 

B(  K1  )u(ad5 


(6.2) 


It  can  be  shown  (c.f.  [35])  that  (6.2)  can  be  solved  by  iteration, 

yielding  the  first  and  higher  order  Born  approximations, 
provided 


X 


2 


r 1 B (r)  1 dr  < 1 . 

0 


(6.3) 


Recent  investigations  in  this  direction  have  been  made  by 
Rorres  ([63])  and  Ahner  ([1]).  V'Je  propose  to  study  (6.1)  by  the 
method  of  integral  operators,  and  shall  show  that  for  X satis- 
fying (6.3)  the  Born  approximations  are  recovered,  whereas  our 
method  also  yields  approximations  in  the  case  when  (6.3)  is  no 

longer  valid.  As  in  the  previous  section  we  assume  that 
2 

B(r)  = 0(e  ) and 

0 < A < 2y(v''2-1)  . (6.4) 


We  restrict  ourselves  to  r 1 and  recall  from  Section  V that 
in  this  case  (c.f.  (5.3)) 


iXz 

e 


+ y b h^^^  (Xr)P  (cos 
In  Bn  ' n ' 
n=0 


6) 


u(x) 


(I+K) 


; r > 1 


(6.5) 
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where  b = b (A)  is  given  by  (c.f.  (5.4)) 
n n 

(I+K)  [(2n+l)i"j^(Xr)+b^(X)h^^^ 

= u^(l)  ^ (1+K)  [{2n+l)i'^j^(Xr)+b^(X)h^^Nxr)  . 


The  quantities  (I+K)e^^^  and  (I+K)h^^Mxr)  can  be  readily 
approximated,  for  small  and  moderately  large  value  of  X,  by 
truncating  the  iterative  process  used  to  construct  the  Icernel 
K(r,s;X)  of  the  operator  K (c.f.  [30]).  Hence  the  problem  is 

to  use  (6.6)  to  approximate  b^(X).  From  (6.6)  we  have  that 
b^(X)  is  an  analytic  function  of  X in  a neighborhood  of  the 
origin  and  has  a Taylor  expansion  of  the  form 

00 

b„(X)  = y b . x’^  . (6.7) 

" k=0 


Hence  a low  frequency  approximation  to  given  by 


bn(X)  ^ b^o^ 


2n+3 


For  larger  values  of  X such  that  (6.4)  is 


still  valid  approximated  by  using  (6.6)  in  con- 

junction with  the  asymptotic  expansions  6f  j^(Xr)  and  h^^^(Xr). 


Returning  now  to  the  low  frequency  approximation  b (X)  % b X 

n nQ 

we  note  that  (c.f.  Section  II) 

2 fl 


2n+3 


G(l,s;X)  = - ^ 


s(l+B(s) )ds  + O(X^) 


(6.8: 


K ( 1 , s ; X ) 


00 

sB ( s) ds  + O ( X^ ) 

1 


and  hence  from  (6.6)  we  can  deduce  that 


.n+1  f..n  , 
1 2 n ! 


2n+l  l(2n) 


2n+2_  , , j 
s B(s)ds 


(6.9) 


This  result  agrees  with  the  first  Born  approximation  (c.f.  [59]) 
in  the  case  when  X satisfies  (6.3). 


We  now  briefly  consider  the  problem  of  low  frequency 
approximations  to  acoustic  scattering  problems  in  a spherically 
stratified  medium  when  an  obstacle  is  present.  For  the  sake  of 
simplicity  we  restrict  ourselves  to  the  case  when  the  obstacle 
is  a "soft"  sphere  of  radius  one  centered  at  the  origin  and 
B(r)  = 0 for  ^ _>  1 . More  general  obstacles  as  well  as  the 

case  when  B(r)  no  longer  has  compact  support  can  also  be 
treated  with  similar  (although  more  complicated)  results  (c.f. 
[28]).  In  our  simple  case  the  relevant  equations  are 


u(x)  = + Ug(x) 


A^u  + X^(l+B(r))u  = 0 in 


(6.10) 


Tim  r - iXu  = 0 


where  h = fx  : jx)  £ 1)  • The  velocity  potential 
scattered  wave  can  be  found  in  the  form 


u (x)  = y b^(X)  (I+K)h^^MXr)P  (cos  01 
s ~ n ~ n n 


(6.11) 
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where  determined  by  the  method  of  separation  of 

variables  ([11],  [25]).  It  turns  out  that  is  an  analytic 

function  of  X in  a neighborhood  of  the  origin  and  has  a Taylor 
expansion  of  the  form 


b (X)  = b , x”^  . (6.12) 

n k=0 


Note  that  in  the  present  case  ^ zero  of  order  2n+l 

at  the  origin,  whereas  in  the  case  of  scattering  without  the 

presence  of  an  obstacle  b (X)  had  a zero  of  order  2n+3 

n 

(c.f.  (6.7)).  In  (6.12)  it  can  be  seen  after  a short  calculation 
([11],  [25])  that  b^Q  is  independent  of  B(r)  and 


ni  n 


B(s)  [s^^'^^t  s”^'^-2s]ds 

1 


(6.13) 


where  y ^ and  y ^ l^nown  constants  (independent  of  B(r)). 

n n 

Hence  for  low  frequencies  the  scattering  due  to  the  obstacle 

dominates  the  scattering  due  to  the  inhomogeneous  medium  and  the 

first  low  frequency  approximation  that  talces  the  inhomogeneous 

medium  into  account  is  b (X)  ^ (b  n+h  ,X). 

n nU  ni 

The  analysis  in  this  section  can  also  be  used  to  investigate 
the  inverse  scattering  problem  of  determining  B(r)  from  a 
Icnowledge  of  the  far  field  pattern  defined  by 
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li^ 


f(0;A)  = lim  re  (x) 

r->oo  s ~ 

(6.14) 

00 

= Y ^ (A)P  (cos  0). 

A n n 


We  restrict  ourselves  to  the  case 

Theorem  6.1;  The  set  |r^^^^+r 

L^[l,a] . 


of  problem 

-2rj  Q is 
J n=0 


(6.10)  . 

complete 


in 


Proof ; [11]  , [25]  . 


If  we  recall  now  the  fact  that 
c(r)  is  the  speed  of  sound  and  c^  ■ 

we  have  from  (6.13)  and  Theorem  6.1 


B(r)  = 


-1  where 


[C  (r) 

= lim  c(r)  is  assuiaad  known 
the  following  result: 


t 


Corollary:  For  problem  (6.10)  the  far  field  pattern  f(0;A), 

O<0£Tf,  0<X<Aq  (where  Aq  is  an  arbitrary  positive 
number)  uniquely  determines  the  speed  of  sound  in  the  inhomogeneous 
medium. 

( 2n+2  ”2n  2 

By  orthonormalizing  the  set  j r +r  '^-2r^  _q  in  L [l,a] 
and  using  (6.13)  to  compute  the  Fourier  coefficients  we  can  also 
obtain  approximations  to  B(r).  However,  since  tends  to 

infinity  as  n tends  to  infinity  this  approximation  procedure 
is  unstable. 


I 


J 


1 
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Open  Problem:  Give  a stable  method  for  constructing  B(r)  from 
f ( e ; A ) . 

Open  Problem:  Consider  the  problems  treated  in  the  last  two 
sections  for  the  case  v;hen  the  source  is  located  a finite  distance 
from  the  scattering  obstacle  instead  of  being  a plane  wave  coming 
in  from  infinity. 

VII.  Heat  Conduction  in  Two  Temperatures. 

In  a paper  on  continuum  thermodynamics  by  Gurtin  and  Williams 
([46])  it  was  shown  that  under  certain  physically  reasonable 
hypotheses  one  could  consider  heat  conduction  being  governed  by 
one  temperature  v;hereas  heat  supply  by  another.  It  was  then 
shown  that  for  an  extremely  general  class  of  simple  materials 
these  two  temperatures  turn  out  to  be  equal.  However,  in  the 
case  of  a non-simple  material,  in  particular  one  in  which  the 
thermodynamic  quantities  depend  on  the  conductive  temperature 
and  its  first  two  spatial  derivatives,  Chen  and  Gurtin  showed  that 
this  was  no  longer  the  case  ([19]).  In  particular  for  an  isotropic 
material  the  linearized  version  of  the  energy  equation  takes  the 
form 


c = kA.j)  I-  ca  A(l>  + q(x,t)  (7.1) 

where  ()j(x,t)  is  the  conductive  temperature,  c the  specific 
heat,  k the  conductivity,  and  a is  the  temperature  descro- 
pency  relating  the  conductive  temperature  to  the  thermodynamic 
temperature  by  the  rtilation 


i 
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0 = (})  - aAtf'  . (7.2) 

From  the  second  law  of  thermodynamics  we  have  the  fact  that 
a ^ 0.  (7.1)  is  a particular  example  of  an  equation  of  pseudo- 

parabolic or  Sobolev  type.  Such  equations  have  been  the  focus  of 
a considerable  amount  of  interest  in  recent  years,  and  the 
reader  is  referred  to  [67]  and  [68]  and  the  references  contained 
in  thses  papers  for  further  information.  In  this  section  we  shall 
restrict  ourselves  to  (7.1)  and  for  the  sake  of  simplicity  set  all 
the  constants  appearing  in  (7.1)  equal  to  unity.  Our  primary 
concern  is  the  solution  of  initial-boundary  value  problems  for 
(7.1)  and  we  first  consider  the  case  of  one  space  variable.  We 
begin  by  introducing  the  idea  of  a Riemann  function.  Let  L and 
M denote  the  operators  defined  by 


respectively.  Then  the  Rieniann  function  v(x,t;C,T)  for 
Ij[u]  =0  is  defined  by 


MlvJ  = 0 
V (C,t;J.,T)  = 

^ (7.4) 

V ( ^ , t ; 4 , r ) = 0 

v(x,T;f,,T)  = 0 


and  can  be  constructed  by  iteration  ([10],  121]).  If  u(x,t)  is 

a solution  of 


L[u]  = q(x,t) 


u ( 0 , t)  = 0 
v.^(0,t)  = git) 
u ( X , 0 ) --  0 


(7.5) 


where  q(x,t)  is  continuous  and  g(t)  continuously  differentiable 


then 


u (x,  t) 


[g'  (t ) (0 , T ;x, t) +g ( t ) (0 , T ; X, t) j dx 


+ q(C,T)v  (5,T;x,t)d5dT  . 

J 0 J 0 ^ 


(7.6) 


Now  suppose  we  want  to  construct  a solution  to  the  initial-boundary 


value  problem 


L[u]  = q(x,t) 

u ( 0 , t)  = 0 
u (x, 0)  =0 
u (a , t)  = 0 


(7.7) 


(We  note  that  the  case  of  non-homogeneous  initial-boundary  data 
can  easily  be  reduced  to  tliis  case)  . From  (7.6)  we  have  the 
following  integral  equation  for  g(t)  = u^(0,t): 

ft 

Y(t)  = g(t)v  (0,t;a,t)  + { v^  ( 0 , x ; a , t ) -v  ( 0 , i ; a , t) ] g ( x ) dx  (7.8) 


where  Y(t)  is  a l<nov;n  function. 


Lemma  7 . 1 : v^ ( 0 , t ; a , t ) / 0 . 

Proof  ( [ ] 0 ) , 1 23  ) ) : From  (7.4)  v;o  have  that  w(x)  = Vj.(x,t:a,t) 

satisfio.s  ii_  -u  U.  lint.  II  (a)  - 0 ^ind  lu^nci.'  ii(0)  = (J  implies 

X X 
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U (x)  = 0 for  0 _<  X £ a which  implies  y'  (x)  = 0 for 
0 £ X £ a,  and  this  is  a contradiction  since  from  (7.4)  y' (a)  = -1. 

We  can  now  conclude  that  the  Volterra  integral  equation 
(7.8)  is  invertible.  Hence  if  we  solve  (7.8)  for  g(t)  and 
substitute  this  back  into  (7.6),  we  have  the  (unique)  solution 
to  (7.7)  : % 

Theorem  7.1;  There  exists  a unique  solution  of  the  initial-boundary 
value  problem  (7.7). 

Proof:  [10],  [21]. 

We  note  in  passing  that  Theorem  7.1  can  easily  be  extended 
to  treat  the  case  of  pseudoparabolic  equations  in  one  space 
variable  with  variable  coefficients  ([10],  [21]).  The  concept  of 
a Riemann  function  for  pseudoparabolic  equations  can  furthermore 
be  extended  to  the  case  of  pseudoparabolic  equations  in  two  space 
variables  ([10],  [22],  [45]). 

We  now  turn  our  attention  to  the  case  wlien  (7.1)  is 
defined  in  a cylindrical  duu\ain  U x (0,T)  where  D n £ 2. 

For  simplicity  we  again  assume  that  all  the  constants  appearing 
in  (7.1)  are  ciual  to  unity  and  set  n = 3.  (Other  values  of  n 
can  be  handled  with  slight  iiiodifi cations ) . We  assume  that  D is 
a bounded,  simply  connected  domain  with  smooth  boujidary  3D, 
and  consider  the  initiai-l^oundary  value  problem 


r 
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+ A^u  = 0 in  D x (0,T) 

u(x,t)  = f(x,t)  on  9D  X [0,T)  (7.9) 

u (x, 0 ) =0  in  D 

where  ffx,t)  is  a continuously  differentiable  function  prescribed 
on  3d  X [0,T].  We  note  that  the  case  of  non-homogeneous  initial 
data  can  be  reduced  to  a problem  of  the  form  (7.9)  by  first 
solving  an  initial  value  problem  by  means  of  the  Fourier  trans- 
form (c.f.[62]).  From  the  maximum  principle  for  pseudoparabolic 
equations  ([71],  [72])  we  have  that  a solution  to  (7.9),  if  it 

exists,  is  unique,  and  hence  our  problem  is  to  derive  a construc- 
tive method  for  obtaining  the  solution  to  (7.9).  We  shall  do  this 
by  the  use  of  a fundamental  solution  and  the  method  of  integral 
equations  (we  note  that  this  approach  can  also  be  used  for 
pseudoparabolic  equations  with  variable  coefficients  ([64]).  We 
define  a fundamental  solution  of  A^u^-u^+A^u  =0  by  the  formula 
([23]) 


l'(x,t;J:,T)  = 


1 

•ni 


R 


U-1  I =C. 


-yR  + ---2  (t-T) 
1-M 


dy 


-R 


I a^ (x; ^)  ( t-T) 
n-1 


n 


(7.10) 


where  R = |x-f,  | and  a^  (x;fj  = 1,  and  look  for  a solution  of 
(7.9)  in  the  form 

3^r 

()(f,  ,t)  (X  , t ; f, , t)  dud  1 (7.11) 

31) 


u(x,t)  = 


2n 
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where  v is  the  unit  norjnal  pointing  into  D and  p(5,t)  is  a 
continuous  density  to  be  determined.  Then  from  the  known  dis- 
continuity properties  of  double  layer  potentials  for  metaharmonic 
functions  we  have  from  (7.11)  that 


(7.12) 


4- 


ao 




9 V a T a t 


(x,t;C,T)dsdT 


where  the  V are  known  constants  (c.f.  [23]).  The  integral 

n 

equation  (7.12)  is  of  the  form  (in  the  Banach  space  C^OD  x [0,T])) 


f^  = (I-T+L.+L-)  p (7.13) 

where  (I-T)~^  exists  (from  the  theory  of  metaharmonic  potentials) 
and  and  are  Volterra  operators.  Hence  in  C^OD  x [0,T]) 

(I-T+Lj^+L^)  exists  and  p = (I-T+L^+L2) 

Opun  Problem;  Consider  oblique  derivative  problems  and  problems 
defined  in  domains  with  moving  boundaries  for  pseudoparabolic 
equations . 

Open  Problem;  Consider  problems  without  initial  conditions  for 


pseudoparabolic  equations  (c.f.  [731  for  the  case  of  the  heat 

equation . ) 
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VIII.  Inverse  1 in  the  Theory  o£  Heat  Conduction . 

In  this  section  we  shnll  consider  two  of  the  classical 
inverse  problems  for  the  heat  equation,  the  problem  of  solving 
the  heat  eq\iation  backv;ards  in  time  (or  the  final  value  problem 
for  the  heat  equation)  and  the  inverse  Stefan  problem.  We  shall 
restrict  ourselves  to  the  case  of  the  heat  equation  in  two  space 
variables.  Other  inverse  problems  for  the  heat  equation  are  also 
common,  and  the  reader  is  referred  to  [2]  and  the  references  at 
the  end  of  [61]  for  further  information.  We  shall  first  consider 
the  case  of  the  final  value  problem  for  the  heat  equation  and 
briefly  indicate  how  an  approximate  solution  can  be  obtained 
through  the  nethod  of  quasireversibility  ([54],  [58])  and  the 

theory  of  pseudoparabolic  equations  as  outlined  in  Section  VII. 
The  final  value  problem  for  the  heat  equation  can  be  formulated 
as  the  problem  of  determining  the  temperature  u(x,y,t)  in 
D x (0,1),  where  D is  a bounded  simply  connected  domain  with 
smooth  boundary  3D,  from  the  equations 


u -ru  = u in 

D (0,1) 

(8.1a) 

xx  yy  t 

u ( X , y , t)  = 0 on 

3D  X (0,1) 

(8.1b) 

u ( X , y , 1 ) = 

f (x,y) 

(8.1c) 

where  f(x,y)  is  a prescribed  function  vanishing  on  3D.  In 
physical  terms  we  are  asking  how  must  tlie  body  D be  heated  in 


order  to  have  a proscribed  temperature  f(x,y)  at  time  t = 1. 
This  problem  i.s  improperly  posed  in  the  sense  that  a solution  does 
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not  in  general  exist,  even  for  f(x,y)  analytic,  and  if  a 
solution  does  exist,  it  does  not  depend  continuously  on  the 
"initial"  data  f(x,y)  (c.f.  [61],  [69]).  However  if  instead 

of  (8.1a)  we  consider  the  pseudoparabolic  equation 

3A2U^+A2U-u^  = 0 (8.2) 

where  3 is  a small  positive  constant,  we  have  from  Section  VII  that 
(8.2),  (8.1b),  (8.ic)  is  a properly  posed  problem  and  can  be 

solved  by  the  method  of  integral  equations  (if  the  Fourier  trans- 
form is  first  used  to  transform  (8.2),  (8.1b),  (8.1c)  to  a non- 

homogeneous  boundary  value  problem  with  homogeneous  initial 
conditions) . The  question,  of  course,  is  what  relation  does  the 
solution  of  (8.2),  (8.1b),  (8.1c)  have  with  the  solution  (if  it 

exists!)  of  (8.1a),  (8.1b),  (8.1c)?  The  answer  to  this  question 

is  contained  in  the  following  theorem  of  Ewing: 

Theorem  8.1:  Let  u(x,y,t)  be  a solution  of  (8.1a),  (8.1b) 

such  that  ||u(x,y,l)  - f(x,y)!|  , < e,  |]u{x,y,0)|l  , < M,  where 

L L 

e,M  are  positive  constants,  and  let  v(x,y,t)  be  a solution  of 
(8.2),  (8.1b),  (B.lc).  Choose  3 = l/log(M/E).  Then  for  every 

t > 0, 

((u-v((  2 - 0 ( (-log  (t/M)  ] ‘^) 


where  the  constant  implicit  in  the  "0"  notation  depends  on  t. 


r 


Proof : [371. 


Note  that  the  constant  M in  Theorem  8.1  can  be  determined 


in  practice  from  a rough  estimate  of  the  maximum  temperature 
attainable  in  the  body  D during  the  time  interval  [0,1]. 


For  other  metliods  for  solving  the  heat  equation  backwards 
in  time  see  [61]  and  [7]. 

We  now  turn  our  attention  to  the  inverse  Stefan  problem 
for  the  heat  equation  in  two  space  variables  ([11],  [20]).  Assume 

that  a bounded  simply  connected  domain  D with  boundary  8D  is 
filled  with  ice  at  0°C.  Beginning  at  time  t = 0 a non-negative 
temperature  y = 'i'(x,t,t)  (where  Y(x,y,0)  = 0)  is  applied  to  3D. 
The  ice  begins  to  melt  and  we  shall  let  the  interphase  boundary 
r(t)  between  ice  and  water  be  described  by 

r(t)  = { (x,y)  : <(’(x,y,t)  = 01  , with  the  water  lying  in  the  region 

<I>(x,y,t)  < 0.  The  ^rft) 


water 


^ / 


-ice  at  0°C 


(normalized)  equations  determining  the  temperature  u(x,y,t)  of 


the  water  are 


where  v is  the  normal,  with  respect  to  the  space  variables,  that 
points  into  the  region  ‘l>{x,y,t)  < 0,  and  the  gradient  is  taken 
with  respect  to  the  space  variables.  The  inverse  Stefan  problem 
is  to  find  y'(x,y,t)  given  '5>(x,y,t),  i.e.  how  must  D be 
heated  so  that  the  ice  nielts  along  a prescribed  path?  This 
problem  is  also  improperly  posed  in  the  sense  that  in  general  a 
solution  does  not  exist,  even  if  >i'(x,y,t)  is  analytic.  In  the 
following  we  want  to  obtain  sufficient  conditions  on  <I)(x,y,t) 
such  that  a solution  docs  exist  to  the  inverse  Stefan  problem. 

Our  approach  also  yields  a constructive  method  for  obtaining 
Y(x,y,t)  {[11],  [20])  but  v;e  shall  not  pursue  this  here. 

For  each  fixed  t,  t e [0,TJ,  let  z = ip(Cft)  map  the 
unit  disc  conformally  onto  a domaiti  such  that  D and 

define  0(f,  ,t)  = (['(C,!)  for  t real.  We  shall  sliow  that  a 
solution  to  the  inverse  Stefan  problem  e.xists  if 


(x  , y , t)  = Im  ^ ( z , t ) 

= yj-  [ !]>  ^ ( Z , t ) -^  ( Z , t ) ] , 


(8.4) 


It  is  assimicd  that  if((S,t)  depends  analytically  on  the  parameter  t 
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Note  that  under  the  assumption  (8.4)  a local  solution  to  (8,3a), 

{8,3c)  exists  by  the  Cauchy-Kowalews)^ i Theorem.  Our  problem 

is  to  establish  the  existence  of  a global  solution,  i.e.  to 

analytically  continue  the  solution  of  a non-characteristic 

* 

Cauchy  problem  for  the  heat  equation.  Let  z = x+iy,  z = x-iy 
^ 

(note  that  z = z if  and  only  if  x and  y are  real)  and 


define 


U ( z , z , t ) 


u 


* 

z+z 

2 


* 


" 2i 


t 


V(  z , z , t) 


1 

t-T 


exp 


(z-U  (z  -C) 
A [t-T ) 


(8.5) 


I 

t 


★ * 

Then  U(z,z  ,t)  and  V(z,z  ,t)  are  solutions  of 


LlUl  5 . 0 


(8.6) 


MIVI  = . 0 


respectively  for  t ^ t.  Let  t be  real  and  for  t on  the  circle 

I t-T I =6,  6 > 0,  let  G(t)  be  a cell  whose  boundary  consists 

-1  - * 

of  a curve  C(t)  ly’ing  on  the  surface  ip  (z,t)  = ip  (z  ,t)  and 

line  segments  lying  on  the  characteristic  planes  z = f,  and 

★ — — 

z = f,  respectively  which  join  the  point  (C(f)  to  C(t). 

A ★ 

Integrating  V[,[Ul-UM(V]  over  the  torus  {(z,z  ,t):(z,z  ) c G(t), 

I t-T I = 6 } gives 


u(r. 


T) 


4 II  i 


★ 

ivtj  ci-’-Vn  *dz  )dt 

Jc(t.)  ='  '■ 

It-, 


) residues 
n=0 


(8.7) 
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where  and  U^*  can  be  computed  from  the  Cauchy  data  (8.3c) 

(c.f.  [11],  [20]).  Due  to  (8.4),  it  can  be  shown  {[11],  [20]) 

that  i)  as  defined  by  (8.7)  is  analytic  in  a region 

containing  D x [0,T]  and  (8.7)  give.s  the  soJution  of  the  inverse 
Stefan  problem  if  we  evaluate  it  on  3D  x [0,T].  Note  that  in 
order  to  get  a physically  meaningful  solution  of  tlie  inverse  Stefan 
problem  we  assume  tliat  ,'(x,y,t)  = 0 on  3D  x [0,T]/^  $ ^ 0 and 
choose  the  conformal  mappings  <[-(^,t)  such  that  u(x,y,t)  ^ 0 
for  4>  < 0 . 

P'or  a partial  extension  of  these  results  to  parabolic 
equations  in  tlirec  space  variables  see  [70]. 

Open  Problem:  Can  the  inverse  Stefan  problem  be  approximated  by 
appropriate  solutions  of  pseudoparaboiic  equations? 

Open  Problem:  What  are  necessary  conditions  for  r(t)  to  be  an 
interphase  boundary  for  tlie  Stefan  problem?  For  a conjecture  in 
the  case  of  one  space  variable  see  [48]  . 

IX.  Hungc's  Theorem  for  Parabolic  Equations  in  Two 

Space  Variables . 

In  this  las),  section  we  shall  outline  liow  the  method  of 
rntegrai  oper.itors  can  be  used  to  obtain  a version  of  Runge's 
Tlicorem  for  parabolic  equations  in  two  space  variables.  Our 
approach  is  a general  i^ation  of  tliat  used  by  Bergman  ([3])  and 
Vekua  ([75])  to  derive  a Uunge  Theorem  for  elliptic  equations  in 
two  independent  variable.s,  ami  the  new  problems  created  in  the 
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present  case  are  that  classical  solutions  of  parabolic  equations 
with  analytic  coefficients  are  not  necessarily  analytic  in  the 
time  variable,  and  in  the  case  of  analytic  solir  ions  we  are 
forced  to  work  v;ith  analytic  functions  of  two  complex  variables 
instead  of  one  complex  variable  as  in  the  case  of  elliptic 
equations.  A further  problem  is  the  need  to  develop  an  analogue 
for  parabolic  equations  of  the  Riemann  function  for  hyperbolic 
equations.  Details  of  the  results  in  this  section  can  be  found 
in  [11]  and  [19].  For  Range's  Theorem  for  parabolic  equations 
in  one  space  variable  see  Theorem  2.1  and  Theorem  3.3  of  these 
lectures . 

We  consider  the  parabolic  equation 

u +u  +a(x,y,t)u  +b(x,y,t)u  +c(x,y,t)u  = d(x,y,t)u  (9.1) 

defined  in  a cylinder  D x (0,T)  where  D is  bounded  simply 
connected  domain  (we  assume  without  loss  of  generality  that  D 
contains  the  origin)  and  we  make  the  assumption  that  the 
coefficients  of  (9.1)  are  real  valued  for  x,y  and  t real  and 
are  entire  funcitons  of  their  independent  complex  variables.  We 
alsp  assume  that  d(x,y,t)  >0  in  D x (0,T).  Let  D be 

simply  connected  and  compact,  6 > 0,  and  D 3 Dq  where 

is  simply  connected  and  such  that  is  analytic. 


bL 


I 


I 


S 
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Theorem  9.1:  For  every  e 
Uj^(x,y,t)  of  (9.1)  in 


> 0 there  exists  an  analytic  solution 
X (0,T)  such  that 


_ max  1 u- -u I < G . 

D^x[5/2,T-6/2] 


Proof  ( [11 ] , [ 19 ] ) : From  the  Weierstrass  approximation  theorem 

and  the  maximum  principle  we  can  approximate  u(x,y,t)  by  a 
solution  Uj^(x,y,t)  assuming  analytic  boundary  data  on  3D. 

The  Theorem  now  follows  from  the  fact  that  solutions  of  (9.1) 
assuming  analytic  Dirichlet  data  on  an  analytic  boundary  are 
analytic  (c.f.  [39]). 


Now  malce  the  change  of  variables 


z = x+iy 
* 

z = x-iy 


(9.2) 


2 

mapping  (t  , the  space  of  two  complex  variables,  onto  itself. 
Then  (9.1)  becomes 


L[U]  = U *^•A(z,  z*  , t)U. +B  (z,  z*  , t)  U * + C ( z , z*  , t ) U-D  ( z , z*  , t ) U = 0 
z z z z z 

(9.3) 

* 11 
where  U(z,z  ,t)  = u(x,y,t),  A = ^(a+ib) , B = ^(a-ib) , 

C = c,  and  D = d. 

•k 

Theorem  9.2:  Let  Uj^(x,y,t)  = U^(z,z  ,t)  be  an  analytic  solution 

of  (9.1)  in  a neighborhood  of  the  origin.  Then  there  exists  a 
* 

)cernel  E(z,z  ,t,T,s)  (whicli  is  independent  of  Uj^(x,y,t))  and 
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an  analytic  function  f(z,t)  such  that 


u^ (x,y , t) 


= Re  P{f} 


= Re  - A(z,o,t)do|  . 


! t-T I =6 


~ ' 2 H T 

E (z  , z , t , T , s) f (z/2 (l-s  ),t) 

-1  ^ 2 

v-  1-S 


Furthermore,  U^(z,0,t)  have  the  same  domain  of 

* 

regularity.  E(z,z  ,t,T,s)  is  an  entire  function  of  its  independent 
complex  variables  except  for  an  essential  singularity  at  t = t. 

Proof  ([11],  [18],  [19]);  Substitute  P{f}  into  L[U]  = 0 and 

integrate  by  parts  to  obtain  a differential  equation  and  initial 

* 

conditions  satisfied  by  E(z,z  ,t,T,s).  This  initial  value 
problem  can  be  solved  by  iteration. 

By  using  the  operator  P we  can  construct  the  complex 
Riemann  function  for  L[U]  = 0 ([11],  [19],  [49])  which  is 

defined  as  the  solution  of  the  singular  initial  value  problem 


M[R]  5 R - AilJI  + CR  + 1^  (DR)  = 0 

2^  o ^ 32 


R (z , z , t ) 


Z 


1 M * 

— cxp^  B(o,^  ,t)da 


(9.4) 


R ( z , z , t ) 


cxp(  A ( f, , o , t)  do  > . 
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Example : For  the  heat  equation  we  have  (see  Section  VIII) 


R(z,z  ,t) 


t-T 


exp 


(2-U 

4 (t-T)  7 


The  complex  Rieiuann  function  is  needed  in  order  to  obtain  the 


following  theorem: 


Theorem  9.3:  Let  (x,y,t)  be  a solution  of  (9.1)  such  that 

Uj^(x,y,t)  is  analytic  for  (x,y)  e D^,  t e 6 wliere  (7  is  a 

★ 

neighborhood  of  [6,T-6].  Then  U-^{z,z  ,t)(=  \i^(x,Y,t))  is 
analytic  in  x x (5  v/here  = { z : z e } , 

* * -k 

= {z  ; z e D^}  . 


Proof  ( [11  ] , [ 19  ] ) : Use  Stokes  Theorem  to  reptos-.nt  u^(x,y,t) 

rn  terms  of  the  complex  Rieraann  function,  v.-hero  the  domain  of 


integration  is  an  appropriate  three  dimensional  torus  situated  in 

3 

(E  , the  space  of  three  complex  variables.  The  regularity  of 
* 

U^(z,z  ,t)  now  follows  from  the  regularity  properties  of  the 
complex  Riemann  function,  which  is  an  entire  function  of  its 


independent  compJex  variables  except  for  an  essentia)  singularity 


at  t = T . 


2 2 “ 1 

Example ; u^(x,y)  = (Itx  +y  ) ''  (this  is  not  a solution  of  (9.1) 

for  any  choice  of  t'ne  coe f L icients  1 ) is  real  analytic  for  all  x 
and  y,  in  particular  in  = {(x,y)  : .x  -ty  <2}  . However  as  a 

■k  k 

function  of  z and  z , U(z,z  ) = u(x,y)  is  singular  in 

* * 

X 0^  , for  oxampJe  at  (i,i)  r Dj^  x . 
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Theorem  9.4  (Range's  Theorem):  Let  u(x,y,t)  be  a solution  of 
(9.1)  in  D X (0,T) . Then  for  every  c > 0 there  exists  an 
entire  solution  u^(x,y,t)  of  (9.1)  such  that  for  ^1^  ^ 

max  ! u-Up I < e . 
n^x[5,T-6] 

Proof  ( [11] , [19 1 ) ; Theorem  9.1  implies  that  u(x,y,t)  can  be 

approximated  on  x [(5/2,T-6/2  1 by  an  analytic  solution 

Uj^(x,y,t)  of  (9.1).  Theoreins  9.2  and  9.3  imply  that 

= Re  P{f}  where  f(2,t)  is  analytic  in  the  product 
domain  x (f.  Since  product  domains  are  Runge  domains  (c.f.[77]) 

f(2,t)  can  be  approximated  by  a polynomial  f^(z,t)  on 
Dq  X (5'^  where  [6,T-6]  C C <?.  Hence  we  can  choose 
U2(x,y,t)  = Re  for  n sufficiently  large. 

Open  Problem:  Derive  Range's  Theorem  for  second  order  parabolic 
equations  in  n _>  3 space  variables  (c.f.  [32]  and  the  references 

cited  there) . 

Open  Problem;  In  Theorem  9.4  can  Dq  x [6,T-6]  be  replaced  by 
D X [0,T]?  Can  the  cylindrical  domain  be  replaced  by  a domain 
with  moving  boundary? 
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